Origin of the "0.25-anomaly" in the nonlinear conductance of a quantum point contact 
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We calculate the non-linear conductance of a quantum point contact using the non-equilibrium 
Greens function technique within the Hartree approximation of spinless electrons. We quantitative 
reproduce the "0.25-anomaly" in the differential conductance (i.e. the lowest plateau at ~ 0.25 — 
0.3 x 2e 2 /h) as well as an upward bending of higher conductance half-integer plateaus seen in the 
experiments, and relate these features to the non-linear screening and pinning effects. 

PACS numbers: 73.23.Ad,73.63.Rt,71.15.Mb,71.70.Gm 



INTRODUCTION 

A quantum point contact (QPC) is a narrow constric- 
tion of a width comparable to the electron wavelength 
defined in a two-dimensional electron gas by means of 
split-gate or etching technique. Due to quantization of 
the transverse motion electrons can propagate only via 
allowed modes and the low-temperature linear-response 
conductance of the QPC shows a step-like dependence 
on a gate voltage [l[ . When a bias voltage V s d is applied 
between the source and drain electrodes the integer steps 
in the differential conductance N x 2e 2 /h are smoothed 
and gradually transformed into the half-integer plateaus 
(N— -|) x 2e 2 /h, where N = 1, 2,3, ... is a number of chan- 
nels available for propagation @, 0, 0, 0, 0, 0, 0, 0, Gil • 
Many features of the linear and nonlinear response of the 
QPC are by now well understood. However, even after 
20 years that have passed since the discovery of the con- 
ductance quantization, some of important aspects of the 
QPC conductance are not resolved yet and are still un- 
der discussions. One of the prominent examples (apart 
from the famous "0.7-anomaly" [llj) is a so-called "0.25- 
feature" in the non-equilibrium differential conductance 
whose origin is under lively current debate [1, 0, [1(3] ■ 

A theory of the non-equilibrium conductance of 
the QPC predicting the above mentioned half-integer 
plateaus was developed by Glazman and Khaetskii [2|. 
The half-integer plateaus have subsequently been ob- 
served and thou ghtf ully studied by a number of groups 
0, 0,0, 0, 0, 0, MUM ■ The theory of Glazman and Khaet- 
skii J2| and a later more refined approach by Frost et 
al. [5j successfully describe the QPC conductance in the 
regime when the differential conductance Gd > 2e 2 /h. 
However, for Gd < 2e 2 /h instead of the expected plateau 
at (0.5) x 2e 2 /h practically all experiments show a plateau 
at (0.2 — 0.3) x 2e 2 /h (sometimes called as a "0.25- 
feature") [1 0, 0, 0, 0, 0, 0, It has been recently 
argued that the "0.25-feature" corresponds to the fully 
spin polarized current even at zero magnetic field [101 ] . 
This conclusion implies far reaching consequences for 



semiconductor spintronics as it opens up exciting pos- 
sibilities to generate spin polarized current simply by ap- 
plying a source-bias voltage to the quantum wire. How- 
ever, alternative explanations of the "0.25-feature" due 
to the self-consistent electrostatics and non-linear screen- 
ing of the lowest spin- degenerate subband have been ad- 
vocated by other groups [1, 0]. In particular, Kothari 
et al. Q demonstrated that the experimental data are 
well-described by the analytical models of Frost et al. H 
with phenomenologically introduced asymmetric voltage 
drop between the source and the drain. 

A detailed understanding of the QPC conductance is 
of the prime importance because the QPC represents the 
cornerstone of mesoscopic physics, and the conductance 
quantization is a fundamental phenomenon of electron 
transport in low-dimensional structures. The contro- 
versy concerning the origin of the "0.25-feature" out- 
lines the need for microscopic modeling based on the 
self-consistent approaches to the electron interaction and 
non-linear screening free from adjustable parameters. It 
should be stressed that previous phenomenological ap- 
proaches 0, 0, 0, 0] , while providing an important insight 
for interpretation of experiments, are not, however, able 
to uncover a microscopic origin of the observed feature. 

In this paper we present a model within the self- 
consistent Hartree approximation that allows us to de- 
scribe the nonlinear screening and evolution of the con- 
ductance plateaus out of equilibrium and thus uncover 
underlying microscopic origin of the observed features in 
the differential conductance. To solve the Schrodinger 
equation we employ a standard non-equilibrium Green's 
function (NEGF) formalism 12, We demonstrate 

that for G < 2e 2 /h the differential conductance ex- 
hibits ~ (0.25 — 0.3) x 2e 2 /h plateau (as opposed to 
the 0.5 x 2eV/i_plateau predicted by the noninteract- 
ing theories [2J, |5j). We also find that in the regime of 
G > 2e 2 /h the nonlinear screening causes the half- integer 
plateaus to bend upward as V s d increases. Note that this 
bending can be clearly seen in all the reported experi- 
ments [1, 0, 0, 0, 0, 0, 0, H3], but, surprisingly enough, 
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its presence passed without comments (except of a brief 
discussion in Q). Our finding therefore strongly indi- 
cates that "0.25-feature" is not spin-related and is caused 
by the non-linear screening and related pinning of spin- 
degenerate electrons in the QPC. 



MODEL 

We consider a QPC defined by split gates in a GaAs 
heterostructure; see Fig. [TJ The Hamiltonian of the 
whole system (the QPC plus the semi-infinite leads) can 
be written in the form H (r) = —^r^ 2 + V e ff(r), where 
r = (x, y), m* = 0.067m e is the GaAs effective mass. 
The effective potential 



lead constriction lead 



V eff (r) = V conf (r) + V H (r) + V blas (r), 



(1) 



is the sum of the electrostatic confinement (includ- 
ing contributions from the top gates, the donor layer, 
and the Schottky barrier), the Hartree and the bias 



potentials (see 
tial is written 



14| for details) 
in a standard form 



Jdv'n(r') 



The Hartree poten- 

i3, m v H (r) = 

where n(r) is 



4Ttc E r J "~ '"^ 1 Ur-r'| A /| r -r'| 2 +4b 2 

the electron density, e r — 12.9 is the dielectric constant of 
GaAs, and the second term describes the mirror charges 
placed at the distance b from the surface, Fig. [TJ The in- 
tegration is performed over the whole device area includ- 
ing the semi-infinite leads; e.g., the Coulomb interaction 
is included both in the leads and in the QPC regions. 

The Fermi energies Ep in the left (L) and right (R) 
leads are shifted by the applied source-drain voltage 
V S ct, Ep — Ep + eV s d, while there is a linear ramp of 
Vbias{ r ) over the device region [l]| (we set Ep = 0). 
For a finite bias the electric current is calculated as 
Ml=fI dE T (E) [fl D (E) - f^ D (E)] , with T(E) 
being the transmission coefficient and f[^(E) is the 
Fermi-Dirac (FD) distribution in the left (right) leads. 
To calculate T(E), the electron density and the local den- 
sity of states (LDOS) we use the standard NEGF method 
0, EH (see Appensix for the details of our calculations) . 
Having calculated the current I we are in position to cal- 
culate the conductance G = I /V s d and the differential 
conductance Gd — dI/dV s d- The latter we compute by 
increasing the bias voltage slightly and calculating the 
derivative dI/dV s d numerically. 

To outline the role of quantum-mechanical effects in 
the electron-electron interaction in the QPC, we also con- 
sider the Thomas-Fermi (TF) approximation solving the 
standard TF equation to find the effective TF potential 
14j and calculating G and Gd for this potential using the 
NEGF. This approximation does not capture quantum- 
mechanical quantization of electron motion and, there- 
fore, utilization of the TF approximation is conceptually 
equivalent to a one-electron noninteracting approach. 
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FIG. 1: (a),(b) Conductance G of the QPC calculated within 
the TF and Hartree approximations for different bias volt- 
ages Vsd- (c),(d) Resonant energy structure (i.e. the LDOS 
integrated over the geometrical area of the QPC constriction) 
in equilibrium, V s d = 0. (e),(f) The differential conductance 
Gd calculated within the TF and Hartree approximations. 
Traces are taken at different gate voltages with 5 mV step 
(the dashed rectangular in (f) contains curves with 2.5 mV 
step);(g) shows the experientially measured Gd adapted from 
[3| • The inset on the top illustrates a geometry of the QPC 
defined in a GaAs heterostructure. A negative voltage is ap- 
plied to the top gates depleting the 2DEG residing on the dis- 
tance b = 60 nm beneath the surface. The widths of the cap, 
donor, and spacer layers are 14, 36, and 10 nm, respectively; 
the donor concentration is 0.64 ■ 10 24 m 3 . The geometrical 
width and length of the constriction are respectively 60 nm 
and 400nm. Temperature T = 0.2 K. 
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RESULTS AND DISCUSSION 

Figure QTa),(b) shows the conductance G of the QPC 
calculated within the TF and Hartree approximations for 
different source-drain voltages V s d- The parameters of 
the QPC are indicated in Fig. [T]and are chosen close to 
those typically used in experiments (Note that we per- 
formed calculations for shorter QPCs which show the 
same behavior). For zero V s d the conductance shows 
well defined quantized plateaus for both TF and Hartree 
approaches. The latter, however, predicts broader tran- 
sition regions between the plateaus. The reason is the 
energy level pinning effect 14|. This is illustrated in Fig. 
Q] (c), (d) that shows the resonant energy structure in- 
side the QPC constriction (i.e. the position of the peak 
in the local density of states (LDOS) integrated over the 
geometrical area of the QPC constriction). In the TF ap- 
proach the resonant levels sweep past Ep in a linear fash- 
ion. In contrast, the Hartree calculations show pinning of 
the energy levels (corresponding to the one-dimensional 
(ID) subbands in the narrowest part of the constriction) 
to Fp within the energy window ±27rfcsT. Within this 
window the FD distribution < f FD < 1, and thus the 
states are only partially filled. This leads to metallic- 
like behavior when electrons can be easily rearranged to 
screen the external electric field. (For influence of the 
pinning effect on equilibrium transport in quantum dots 
see [lj]; see also [16| for the experimental studies of the 
energy level pinning in the QPC in the magnetic field). 

Out of equilibrium, the energy window eV s d provid- 
ing current carrying states increases as the source-drain 
voltage grows, and the conductance plateaus become 
smeared, see Fig. [I] (a), (b). The plateaus in the conduc- 
tance G completely disappears when eV s d exceeds the ID 
subband energy separation inside the QPC constriction. 
At the same time, half-integer plateaus (N — \) x 2e 2 /h 
appear in the the differential conductance Gd, see Fig. 
Q] (e), (f)). A comparison of the TF and Hartree re- 
sults show two profound differences between the calcu- 
lated Gd- First, the lowest Hartree plateau N = 1 
rs 0.25 — 0.3 • 2e 2 /h as opposed to the 
0.5 x 2e 2 /h plateau predicted by the TF cal- 
culations. Second, all TF plateaus are flat and rather 
independent of V s d, whereas all higher Hartree plateaus 
N > 2 are bent upward as V s d increases. Note that these 
two features of the calculated are clearly seen in all 
reported experiments % i, H, [EH B @, G3 (see Fig. Q] 
(g) for a representative example). 

In order to understand the origin of the above fea- 
tures of the QPC nonlinear conductance let us inspect 
the LDOS inside the QPC region. Let us first concen- 
trate at the first plateau in the differential conductance. 
Figures [2] (c), (e), (g) show the evolution of LDOS as 
V s d is increased calculated within the noninteracting TF 
approach. The enhanced LDOS in the constriction corre- 
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FIG. 2: The charge density and the LDOS of the QPC calcu- 
lated in the Thomas-Fermi and Hartree approximations (left 
and right panels respectively) for the first half-integer plateau 
for different V s d- The corresponding gate voltages V g are 
marked by arrows in Figs. [lja),(b),(e),(f). The effective po- 
tential Eq. ([l]) is plotted by the dashed lines. Solid slanted 
lines denote the bias potential profile Vbias- Insets show the 
current profiles, T(E) [f[ D (E) - f£ D (E)]. 



sponds to the position of the bottom of the lowest prop- 
agating subband. In the TF approximation the effective 
confinement potential is symmetrically distributed rela- 
tive to Vbias (that ramps linearly along the device). Be- 
cause of this the ID subband touches Vbias at the QPC 
center (at the energy E = Ep — ). As a result, the 
electrons injected from the left lead in the upper half of 
the eV sd window (Ep < E < E- — pass through 

the QPC with the unitary probability. However, the elec- 
trons in the lower half of the eV s d window experience a 
potential barrier and hence are reflected back (see partial 
current profiles, T(E) [f[ D {E) - f£ D {E)] , in small in- 
sets in Figs. [2](c), (e), (g)). Thus, the electrons injected 
from the left lead give rise to the conductance of the half 
of the conductance unit, Gd = 0.5 x 2e 2 /h. For electrons 
moving in the opposite direction, from the drain to the 
source electrode, there is no available channel to propa- 
gate and all of them are reflected. 
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FIG. 3: TF and Hartree LDOS in the QPC for different 
V s d- (The same as in Fig. [2] but for the second half- integer 
plateau in the differential conductance, see arrows at Figs. 
[TXa),(b),(e),(f).) 



A character of electron transport changes dramatically 
when interaction is included at the quantum-mechanical 
level. Figures [U (d), (f), (h) show the LDOS inside the 
constriction calculated within the Hartree approximation 
for the first half-integer plateau where Gd ~ 0.3 x 2e 2 /h. 
With one partially propagating mode the electron den- 
sity inside the constriction is low and the screening is 
rather week, and hence the electron interaction strongly 
modifies the potential profile in comparison to the sym- 
metric TF distribution. The Coulomb charging pushes 
up the upper ID subband inside the QPC constriction 
to the top of the eV s d window near the source con- 
tact. (It is interesting to note that the LDOS inside the 
QPC out of equilibrium resembles a corresponding self- 
consistent LDOS profile of a resonant-tunneling diode 
[ftI]). Thus, the ID subbands becomes pinned to the top 
of the eV s d window and therefore only a relatively nar- 
row energy interval can supply current-carrying states 
that pass through the QPC (see current profiles in the 
insets to Fig. H (d),(f),(h)). Hence, the QPC con- 
ductance, Gd ~ 0.3 x 2e 2 /h, becomes smaller than a 
half of the conductance unit Go = 2e 2 /h. Our calcula- 
tions provide therefore a microscopic foundation of the 
phenomenological approaches that describe the "0.25- 
anomaly" assuming a nonsymmetric voltage drop inside 
the constriction [§]. 

Let us now turn to higher half- integer plateaus. In this 
case there is at least one propagating state inside the 
constriction, which, in turn, leads to enhanced screen- 
ing. Indeed, despite of the voltage drop between the left 
and the right leads, the Hartree effective potentials V e // 
and the LDOS are practically flat inside the QPC, see 
Figs. [3] (b), (d). This is in contrast to the correspond- 



ing TF results which do not account for screening and 
thus follow the linear drop of Vbi as , see Figs. [3] (a), (c). 
Because of the enhanced screening, at the center of the 
QPC the Hartree ID subbands are situated lower than 
the corresponding TF subbands (i.e below Ep — ), 
and hence the energy window providing the transmitted 
states through the QPC exceeds the half of the available 
energy interval eV B d (see the current profiles in the inset 
of Fig. [3] (b), (d)). As a result, the QPC conductance 
corresponding to the highest subband is larger than a 
half of the conductance unit Go- (Note that all lower 
subbands are fully occupied and thus contribute to one 
conductance unit each). Thus, the enhanced screening, 
which becomes more pronounced as V s d increases, is the 
reason for the upward bending of the higher half-integer 
plateaus. 

Finally, we stress that we utilized a model of spinless 
electrons in the Hartree approximation. The present ap- 
proach can be easily extended to account for the spin ef- 
fect within the framework of the spin-density functional 
theory (SDFT). However, some previous studies ques- 
tioned the reliability of the SDFT for the system at hand 
because of the self-interaction errors of the local spin 
density approximation [3]. Hence, the definite answer 
about the role of the spin in the non-linear conductance 
of the QPC might require approaches that go beyond 
the mean field method used in the present study (e.g. 
quantum Monte-Carlo, etc.). At the same time, an ex- 
cellent quantitative agreement of our calculations with 
the experimental results outlines the dominant role of the 
self-consistent electrostatics and the nonlinear screening 
and strongly indicates that the "0.25-feature" is not spin- 
related. 

To conclude, using NEGF formalism within the 
Hartree model of spinless electrons we reproduced quan- 
titatively the observed features of the nonlinear QPC 
conductance and provide microscopic interpretation of 
the "0.25-anomaly" as well as the upward bending of the 
higher half-integer plateaus in terms of non-linear screen- 
ing and pinning effect. 

This work has been supported by the Swedish Research 
Council (VE). 



APPENDIX: NON-EQUILIBRIUM GREENS 

FUNCTION (NEGF) TECHNIQUE FOR 
CALCULATION OF THE TRANSMISSION 
COEFFICIENT OF THE QPC 

The central quantity in the NEGF is the lesser Green's 
function, G < [l2| . To calculate it one has to find first the 
retarded Green's function, G r , 



(E-H(r))G r (r,r',E) = l, 



(2) 



where E is an electron energy and 1 is the unitary op- 
erator. This equation can be reformulated using the so- 
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lm(E) 



Path 1 




FIG. 4: An integration path used in Eq. 0. Path 2 appears 
when the bias voltage V s< j applied. 



called retarded self-energies of the leads, E^ and E£, 

(E-H - E R (E) - YT R (E)) G r (E) = 1, (3) 

where Hq is the Hamilton operator for the isolated scat- 
tering region (i.e. excluding the leads). (For the sake 
of shortness we will not write explicitly a coordinate de- 
pendence of G r ). ^n(L) are functions with non-zero val- 
ues only at the boundaries with the semi-infinitive leads. 
Coupling the scattering region with leads is described by 
the functions 



iT R (E) = l, r 
iT L {E) = 1 



(E)-X R (E) = 2ilm[X R (E)}, (4a) 
l{E)-Vl{E) = 2ilm[W L {E)]. (4b) 



2 in Fig. [4jb) where f R R\ ^ 1)j ^ s c l° se to the real axis 
and there Eq. ([5]) is used. Along the Path 1 only several 
integration points are needed because the rapid varia- 
tions of Gf q {E) are smeared out when the integration 
path is far from the real axis. This is specially useful for 
the bound states, which give rise to sharp peaks near the 
real axis. On the straight path along the real axis, one 
needs much more integration points and for large source- 
drain voltage it becomes the most time consuming part 
of computation. 

Equations ©-(O are solved self-consistently in an it- 
erative way until a converged solution for the electron 
density and potential (and hence for the total Green's 
function) is obtained. Having calculated the total self- 
consistent Greens functions, the transmission coefficient 
is calculated as 12 1 



T(E) = Tr [T L (E)G r (E)T R (E)G a (E)} 



(8) 



To speed up computation we employ the hybrid recursive 
method working with sin-Fourier transformed Greens 
functions and use the second Broyden method for 
the iterative algorithm 21 1. 



The lesser Green's function in the scattering region is 
related to electron flow from right and left reservoirs and 
is written as 

G<(E) = - if^ D {E) G r (E)T R (E)G a (E) 

- if[ D (E)G r (E)T L (E)G a (E), (5) 

where f R R) are the Fermi-Dirac functions in the right 
(left) lead. This equation has to be used in non- 
equilibrium situations when V s d ^ and f R D ^ /™, 
In equilibrium, when the Fermi functions in both leads 
are identical, Eq. (JSJ) reduces to 



G<(E) = 2f™ L) (E) G r (E) 



(6) 



It is also valid under a bias voltage at energies E for which 
f R D = f[ D (in practice, f R R\ = 1 for those energies). 

In order to calculate the electron density we integrate 
over the electron energy E 



i(r) = -i- / dE Im [G< (r, r , E) 



(7) 



We use both Eqs. (|5|) and ([6]) to perform this integra- 
tion. Gf q (E), Eq. ©, is analytic in the upper half of 
the imaginary plane whereas G < (E), Eq. (JSJ), has poles 
below and above the real S-axis. Thus, for the energies 
when f R R\ = 1 we can use Gf q (E) for which we can 
transform the integration path from the real axis to the 
complex plane [lj| , see Fig. [H where Gf q (E) is a smooth 
function of energy. The rest of the integration (i.e. Path 
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